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X-ray absorption and occultation in LS 5039 
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ABSTRACT 

Gamma-ray binaries are systems containing a massive star and a compact object 
that have been detected up to TeV energies. The high energy emission could result 
from particle acceleration in the region where the stellar wind from the massive star 
interacts with the relativistic wind from a young pulsar. LS 5039 has the most compact 
orbit amongst gamma-ray binaries and its X-ray lightcurve shows a stable modulation 
synchronized with the orbital period. Photoelectric absorption of X-rays in the O star 
wind and occultation of the X-ray emitting region by the massive star can alter the 
X-ray lightcurve and spectrum along the orbit. Yet, the X-ray spectrum and lightcurve 
of LS 5039 do not show intrinsic absorption or X-ray eclipses. We study these effects 
in the framework of the pulsar wind scenario as a function of the binary inclination 
angle, the stellar wind mass-loss rate and the size of the X-ray emitter. An extended 
X-ray emission region > 3i?+ appears necessary to reconcile the pulsar wind scenario 
with observations. 

Key words: binaries: close - binaries: eclipsing - stars: individual: LS 5039 - stars: 
mass-loss - gamma-rays: stars - X-rays: binaries 



1 INTRODUCTION 

Gamma-ray binaries are systems containing a massive star 
and a compact object that emit most of their power at en- 
ergies above 100 MeV. They have been detected up to very 
high energy (VHE) gamma-rays and are thus sites of particle 
acceleration up to multi-TeV energies. Ther e are three estab- 
lished gamma-ray binaries PSR B1259-63 (Aharonian et al. 

1 ° * — ■■ ■ i !^ 



lished gamma-ray binar ies ran aLZo'd-bd | Aharoma n et al.l 
l2005t ). LS 5039 (lAharonian et al.ll200rjl : lAbdo et al.ll2009b ) 
and LS 1+61 303 (|Albert et alj 1200a Ubdo et al.l l2009ah 
and a candidate binary HESS J0632+057 IjAharonian et al.l 
l2007h . 

The VHE emission in gamma-ray binaries is thought 
to be due to Compton upscattering of UV photons 
from the massive star by energetic electrons. Electron 
acceleration could take pl ace either in a relativistic jet ( mi- 
crogl ias ar scenario. iRomero. Christiansen fc Orellanal 



20051: iParedes Bosch-Ramon fc Romer J ~ l2006t 

Dermer fc Bottched l2006f > or in the shocked wind of a 
young pulsar (pu lsar wind scenario, iMaraschi fc Treves! 
ll98ll : lDubu"sll2006h . where the shock results from the inter- 
action between the pulsar wind and the stellar wind of the 
companion star. The latter model is known to be operating 
in case of PSR B1259-63 llTavani. Arons fc Kaspil 1 1994 
iKirk. Ball fc Skiaeraasenl Il999h . The nature of the compact 
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object and the electron acceleration site are uncertain in 
the remaining gamma-ray binaries. 

Gamma-ray binaries have also been observed in X-rays 
where their properties differ from those of high-mass X-ray 
binaries (HMXB). The X-ray spectrum of gamma-ray bi- 
naries are power-laws but show no apparent cutoffs up to 
hundreds of keV. The X-rays are thought to be due to non- 
thermal synchrotron or inverse Compton emission. Gamma- 
ray binaries also do not display X-ray outbursts and state 
transitions as usually seen in accreting binaries. 

LS 5039 is the most compact gamma-ray binary, com- 
posed of an unknown compact object in a 3.9 day orbit 
around a 06.5V star l|Casares et al.ll2005h . It has a regu- 
lar behavior in both gamma-rays and X-rays. The X-ray 
lightcurve shows an orbital modulation w i th a remarkable 
long- t erm stability dKishishita et al.l 120091 : iTakahash i et al.1 
2009; Hoff mann et al.l 12009] ). The modulation is present at 
low (1-10 keV) and medium (10-40 keV) X-ray energies, 
with minimum and maximum flux at superior and inferior 
conjunction respectively. The location of the minima and 
maxima suggest the modulation is a geometrical effect re- 
lated to the orientation of the binary with respect to the 
observer, rather than due to physical changes taking place 
in the shocked winds region when the compact object travels 
on its elliptical orbit. 

X-ray absorption in the stellar wind and occultation of 
the X-ray emitting region by the massive star both result 
in orbital modulations with the correct phases for flux min- 
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Figure 1. The schematic 2D illustration of the model where the 
stellar wind and pulsar wind collide. The sizes of the pulsar and 
star are not to scale. 



imum and maximum. However, for LS 5039 there is no evi- 
dence for an absorption excess due to the st ellar wind at any 
orbit al phase, even at superior conjunction ((Takahas hi et all 
2009). This has been used to argue that the X-ray emit- 
ting region must be far out from the system, where the 
column density of materia l crossed by the X-ra ys is small 
i|Bosch-Ramon et alj [20071 ; iBosch-Ramonl |201Ch . Here, we 
re-examine this question in detail, taking into account the 
3D geometry of the interaction in the pulsar wind scenario 
to compute model lightcurves. This is used to derive con- 
straints on the binary parameters from the lack of absorp- 
tion and occultation signatures in the observed lightcurve. 

The outline of this article is as follows. The model is 
described in fJ3] and applied in 33] to estimate the size of 
the X-ray emitting region based on the observed hydrogen 
column density. The occultation of the X-ray emitting re- 
gion by the massive star is studied in section Finally, Sj5] 
discusses the findings and presents the conclusions. 



2 THE GEOMETRIC MODEL 
2.1 Shape of the shock 

In the pulsar wind scenario, X-rays are emitted by parti- 
cles accelerated in the shock region where the pulsar wind 
and the stellar wind of the massive star collide (Fig.[T]). The 
interaction region is bounded by two termination shocks S 1 * 
and S p . Downstream (zones 2 and 3), two shocked winds are 
separated by a tangential contact discontinuity (CD). Up- 
stream (zones 1 and 4), the winds behave as in the case of 
single star/pulsar. In this example, the stellar wind momen- 
tum dominates over the pulsar wind. The stagnation point 
Rs between the two winds along the line joining the stars is 
found by equating the ram pressures from the two winds: 
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Mv(R s ) = 
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where a is the binary orbital separation, M is the stellar 
wind mass loss rate and E is the pulsar spindown power. The 
stellar wind velocity v at the radial distance r from the star's 
centr e is given by a /3-velocity law jCastor. Abbott fc Kleinl 



v(r) 



ro 



(2) 



where ro = 1 — (vo/Voo) with no the initial wind velocity, 
Vaa is the wind terminal velocity, R+ is the stellar radius and 
f3 ~ 1 is a parameter describing wind acceleration. 

The detailed structure of the shock will depend on the 
momenta of the winds and on radiative cooling of the gas. 
Efficient cooling will tend to collapse the termination shocks 
onto the CD. The shock structure can also be affected by 
mixi ng instabilities at the interface a nd by orbital motion 
(e.g. IStevens. Blondin fc Pollocklll992l ). There is no general 
semi-analytic description of the shock structure but a de- 
scription of the CD can be obtained. The CD is defined as 
the surface where the perpendicular components of the ram 
pressures balance each other p t ± — p p ±. For each point in 
space one can then define a dimensionless parameter 



V{ri) = 



E 



r\ sin 2 i 



Mv(n) 



(3) 



where 9\ an 62 are angles between the line tangential to the 
CD at the given point and the direction towards the star 
or pulsar respectively as illustrated in Fig. [1] 7*1 and ri are 
distances from the star and pulsar respectively. The value of 
77 at the stagnation point (d\ = 62 — t/2) equals 



7/0 = r)(R s ) = 



= 0.05 



R B f 



Mv{R s 
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10 36 erg s - 
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10- 7 M yr- 
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10 8 cm s -1 



The 770 depends on orbital phase via v(R s ) and parametrizes 
the shape of the shock, which is then obtained by solving the 
differ ential equation (Eq. 6 in lAntokhin. Owocki fc Brownl 
120041 ) 



dx 
dy 



arii^yh/v [ri(x,y)] 
ri( x >y)Vv [ri(x,y)} + r%(x,y) 



(5) 



with the initial condition x = R s and y = 0. Numerical simu- 
lations show that this holds reasonably well in the case of the 
inter action with a relativistic pulsar wind l|Bogovalov et all 
120081 ). The solution is a 2D profile of the CD x = f(y). The 
CD is symmetrical with respect to a;-axis, thus it is best 
represented in cylindrical coordinates (p,ip,x) as x = f{p), 
where p — y 2 + z 2 and ip = arctan(z/y). 

An example 3D representation of the model is shown 
in Fig. [2] In cases where 770 < 1 the stellar wind dominates 
over the pulsar wind. When 770 <C 1, the opening angle of the 
CD a < 45° (where a is an angle between x-axis and a line 
tangential to the CD at a large distance from the stagnation 
point) and the CD wraps around the pulsar creating a tail- 
like structure. For 770 > 1 the opening angle a > 90° and 
the CD curves around the massive star. Note that there is a 
maximum value of 770 associated with a minimum realizable 
distance _R™ m between the star and CD. This is approxi- 
mately the place where p+(r) has its maximum. Stable bal- 
ance is lost if the pulsar wind moves beyond this point: then 
the pulsar wind overwhelms the stellar wind and is stopped 
at the star surface. 

The shape of the CD will be correct up to the point 
where the Coriolis force curves the shock structure. Assum- 
ing that the stellar wind is collimating the pulsar wind, this 
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Figure 2. 3D illustration of LS 5039 in the framework of the 
pulsar wind scenario. The plot is in units of orbital separation and 
the massive star is plotted to scale. The dark sphere is the massive 
star, the light grey (green) points represent the CD surface, the 
dark grey (red) points are the points that emit X-rays, whereas 
black points are in the shadow of the massive star. The arrow 
represents the direction to the observer at inclination 60° and at 
the phase of superior conjunction. 



happens at a distance ~ Voo-Porb, or about 5.4 a.u. (130 R*) 
for LS 5039 with the parameters given in £12.41 Closer in, 
the orientation of the surface can also be altered by orbital 
motion: the skew angle s is given by tan s = v OI b/v(a) where 
Worb is the pulsar orbital velocity. 



2.2 X-ray emission 

In the pulsar wind scenario, the X-ray emission is primarily 
due to synchrotron emission in the shocked pulsar wind re- 
gion (zone 3 in Fig. [TJ which is tenuous and optically thin 
to X-rays. The stellar wind region may also be a source 
of X-rays, either because of the kinetic energy dissipated 
in the colliding winds (zone 2) or du e to small-scale insta- 
biliti es in the stellar wind (zone 1) l|Puls. Vink fc Naiarrol 
2008). In both cases the emission is softer (kT ~ 0.5 keV) 
and weaker (10 30 -10 33 erg s" 1 , e.g. lStevens et al.lll992T ) than 
the observed X-ray em ission from LS 5039 (~ 10 33 erg s _1 , 
iTakahashi et al.ll2009T l. Detecting X-ray line emission would 
provide valuable diagnostics of the stellar wind but this 
emission is most likely swamped by synchrotron emission. 

T he cooling paramet er \ = icooi/tdyn of the O star 
wind (|Stevens et al.ll 19921 ) is « 1 in LS 5039, at the limit of 
efficient radiative cooling. Isothermal (radiatively efficient) 
shocks have small widths. The pressure in the shocked pul- 
sar wind will mostly be due to the lowest-energy electrons if 
the particles have a steep distribution. These electrons have 
long cooling times and the shock may be more closely ap- 
proximated as adiabatic rather than isothermal. Complex 
numerical simulations are required to obtain the detailed 
structure and emissivity of the shock region. At this stage, 
we assume that the shocked regions 2 and 3 have a negligible 
width and that the emissivity is uniform over the 2D shock 
surface, which should suffice to capture the geometrical ef- 
fects that we want to investigate. 

The 2D shock surface is approximated by the CD sur- 
face for p £ (0, /5 max ) and zero otherwise, with p ma x the 



cylindrical radius at which a sphere centered at (R S ,Q, 0) of 
radius i? ou t intersects the CD (see Fig. [1} . The size of the 
emitting region is conveniently parametrized by R ou t in the 
following. The total flux from the emitting volume is 



F{E) = [ j E dV = f X j E dS, 
Jv J s 



(6) 



where j'e is the unit volume emissivity, and Aj'e is the con- 
stant surface emissivity. The unit surface element of CD 
equals 



dS 



1 + 



p dp dtp. 



(7) 



2.3 Absorption and occultation 

The uniform emission of X-rays is affected in two, orbital 
phase-dependent, ways. First, X-rays can undergo photo- 
electric absorption as they cross the dense stellar wind. The 
wind density at the distance r from the stellar center is 



n(r) 



M 



(8) 



4:TT(imnr 2 v(r) ' 

where tjih is hydrogen mass while fi = 1.3 is the mean molec- 
ular weight. Second, at some phases and for some inclina- 
tions, parts of the emitting surface of the CD are occulted 
by the massive star and are thus invisible to the observer 
(black regions in Fig. [2} ■ 

In order to estimate the impact of absorption and occul- 
tation on the observed flux, a random sample of TV points 
is uniformly distributed on the surface of CD out to R ou t 
(see Fig. [2] and Appendix A) . Each point corresponds to a 
7V-th part of the emitting surface and has a flux equal 1/N 
(the total flux is normalized to unity) . The observed flux at 
given energy E, inclination angle i and orbital phase <j> is 



1 N 

F(E,i,<j>) = -£exp[-a(E)JV Hj (i,<«] CiM), 



(9) 



3 = 1 



where Nuj is the stellar wind's hydrogen column density ob- 
tained from integration of n(r) along the line of sight from 
point j, a is the photoelectric cross section of the plasma 
which is a function of photon energy E. If the plasma is ion- 
ized, a also depends on plasma ionization parameter which 
varies along the line of sight, is the occultation function 
which value is whenever a line of sight crosses the interior 
of the massive star and 1 otherwise. 

In calculations of TVh the parts of the line of sight 
that cross the pulsar side of the CD (zones 3 and 4 in 
Fig. [TJ are considered to be empty since neither unshocked 
nor shocked pulsar wind have sufficient densities to ab- 
sorb X-rays. For simplicity, we also assume that zone 1+2 
are described by Eq. [5] all the way from the stellar sur- 
face out to the CD i.e. we ignore the density enhance- 
ment in the shocked stellar wind (zone 2). Again ( cj2.2[) . 
detailed numerical simulations would be required to model 
this appropriately. An estimate ca n be derived using the 
thin s hell colliding wind model of ICanto. Raga fc Wilkinl 
( 1996). The surface density of the shock is « 3M y/rj/ (167ra«) 
in the orbital plane, where it is maximum. This corre- 
sponds to a column density across the shock TVh, shock ~ 
7 W 19 M 10 -r M(3 /yrVoma-o.i au w M00 km/s cm " 2 > compared to 
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the wind column density Nh ~ 2 10 21 cm -2 from R+ to in- 
finity (for the same wind parameters). A^h, shock is also about 
30 times less than the inferred wind column density at su- 
perior conjunction (see Sj3}. We conclude that taking into 
account the density enhancement in the shocked region is 
unlikely to change our results significantly. 

2.4 Parameters for LS 5039 

We adopt the following binary parameters for LS 5039: 
^oo = 2.4 x 10 s cm s~\ i?* = 9.3-Rq, M* = 23M and = 
3.9 x 10 4 K. For the wind of the massive star, we adopt /3 = 1 
and u = 2 x 10 6 cm s~ 1 (Eq. O . The value of M in LS 5039 is 
not well constrained and varies from ~ 3 x 10 -8 to 7.5 x 10~ 7 
Mq yr _1 , with values of several 10~ 7 Mp yr -1 fa v oured 
in the most recent litterature (Ku dritzki fc Pulsl |2000|; 



McSwain fc Giesll2002l; IMcSwain et al.ll200j ; ICasares et"aH 
20051 ; ISzalai. Kiss fc Sartvll2010l ). These mass loss rates are 



derived from Ha line fitting, a diagnostic whi ch is known to 
be affected by wind clumping (|Puls et al . 2008). The true M 
could be a factor 2-3 lower than the M estimated from Ha. 
The large range of r\ that we explore covers this uncertainty. 

The binary inclination angle is i > 40° if the binary con- 
tains a neutron star. The orbit eccentricity is e = 0.33, while 
the a ngle of the line of nodes is cj = 236° l|Aragona et al.l 
2009). For these orbital parameters the periastron, apas- 
tron, pulsar superior and inferior conjunction correspond to 
orbital phases cj> — 0, 0.5, 0.045 and 0.67 respectively. The 
respective binary separations are 1.4 x 10 12 , 2.82 x 10 12 , 
1.46 x 10 12 , and 2.59 x 10 12 cm. The skew angle is s = 23° 
at periastron and 8° in apastron. This has a negligible im- 
pact on the conclusions and is not taken into account. 

A strict lower limit on the value of E in LS 5039 is 
~ 10 35 erg s _1 which implies 100% radiative efficiency to 
gamma rays (at a distance of 2-3 kpc). At the other end 
of the scale, the pulsar wind will impact the star surface 
beyond a certain £ max set by M and the binary separa- 
tion at periastron. Such a situation is thought to occur in 
the black widow pulsars where the wind from the lo w-mass 
stellar companion is very weak (|Phinnev et al.|[l988l ). If the 
pulsar wind is dominated by the kinetic energy of the cold e + 
e~ pairs then radiative braking of the pulsar wind by inverse 
Compton losses m ight still allow a shock to form just above 
the O star surface (jCerutti. Dubus fc Henrl2008l ). The mag- 
netic field of the O star may also be suffi cient to hold off 
the pulsar wind (|Harding fc Gaisserlll990T l. However, there 
is no observational evidence that the stellar wind in LS 5039 
collapses due to the relativistic wind momentum impi nging 
on its pulsar-facing hemisphere. IMcSwain et al.l |2004l ) find 
no differences in UV spectra taken at two different orbital 
phas es (0.40 and 0.63 based on ephemeris of lAragona et all 
2009), the latter one close to inferior conjunction of the com- 
pact object where the effect of the O star wind quenching 
by the pulsar wind should be most visible. Assuming the O 
star wind does not collapse implies a maximum 77 max ~ 0.6, 
reached at periastron, for which the CD opening angle is 
~ 62°. For the lower and upper limit on M, the correspond 
ing iJmax varies between 



1.5 x 10 36 erg s" 1 and 3.8 x 10 37 



the pulsar wind. At this value of 770 the CD opening angle 
is only ~ 10°. 

Ther e is also no evidence fo r X-ray ionization of the stel- 
lar wind (Mc Swain et al.l|2004l ). suggesting that the impact 
of ionization on X-ray absorption can be neglected. We con- 
firmed the validity of this assu mption for LS 5039 with the 
photoionization code XSTAR (Kallman 2005). In the simu- 
lation, an optically thin plasma cloud was illuminated by a 
source with a spectrum composed of a power law with spec- 
tral index 1.5 and 1-10 keV luminosity 6.03 x 10 33 erg s" 1 
(|Takahashi et al.ll2009l assuming a 2.5 kpc distance to LS 
5039 from Earth) and a black body (stellar continuum) with 
a temperature T*, radius i?* and luminosity L* = 6.9 x 10 38 
erg s _1 . The measured opacities in the 1-10 keV range 
did not differ from thos e of th e cold plasma described by 
iMorrison fc McCammor] (| 19831 ). The intense UV radiation 
of the massive star does not influence the plasma opacity at 
energies above 0.5 keV. 



3 X-RAY ABSORPTION IN LS 5039 

Observations show no signatures of X-ray absorption by the 
stellar wind in LS 5039. The XMM-Newton and Suzaku spec- 
tra of LS 5039 are absorbed by an equivalent hydrogen col- 
umn density consistent with the Galactic value. The upper 
limit on the intrinsic column density at the phase of superior 
conjunction (where the line of sight column density related 
to the stellar wind is highest) is as low as N„ ax = 2.6 x 10 21 
cm -2 (|Bosch-Ramon et al.ll2007l ; iTakahashi et alj|200gh . 

The photoelectric absorption depends on four param- 
eters, M, i, Rout and r/o- This can be used to estimate 
the maximum allowed Af max for which the line of sight 
column density does not exceed 7V^ ax . For example, for 
a point X-ray source located at the pulsar position, at 
the phase of superior conjunction, and at inclination angle 
i — 40°, the upper limit on the stellar wind mass loss rate 
is Mmax = 7.3 x 10~ 8 Mq yr" 1 . For a point source located 
at the stagnation point, the Af max is lower, since the line of 
sight probes regions of the stellar wind with higher density. 
For an X-ray emitting point source located at R s and for 
which 770 = J?max, then Af max = 2.3 x 10~ 8 M Q yr _1 . 

These values of M max are at the lower end of the M 
scale for the Galacti c O stars, (~ 10~ 8 - 10~ 6 M yr" 1 , e.g. 
iMokiem et~alll2007l ). Af max could be higher if the wind was 
highly ionized but this is unlikely (see ij2.4p . A/ max could also 
be higher if the X-ray emitter is away from the massive star, 
e.g., in the form of a jet, emitting X-rays at a distance > 10 12 
cm w here the stellar wind is rarified (Bosch- R amon et al.l 
2007). Another possibility is that the X-ray emitting region 
is not point-like but large compared to the size of the massive 
star. 

For the measured value of A^g lax and as a function of 
Af max , we estimate the minimum size of the X-ray emitting 
region for which intrinsic absorption would not manifest it- 
self in the X-ray spectrum. The upper limit N™* is based 
on X-ray data fitting which assumed that an entire (point 
or extended) X-ray source is covered by a uniform absorbing 
gas cloud so that 



(4.6 x 10 erg s~ ). In the following calculations we often 
use the value of rjo ~ 0.004 to examine the model predic- 
tions in the case of extreme stellar wind domination over 



F(E) =exp[-tr(£)A%], 



(10) 



where the intrinsic flux equals 1. In our case however, where 
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Figure 3. The minimum size of the X-ray 1 keV emission region 
at superior conjunction. The dots mark the results of calculations, 
the lines are linear fits to the points. The dashed lines correspond 
to r]o = 0.004 at superior conjunction, the solid lines correspond 
to maximum value of 770 = t)mai = 0.6 at periastron. The grey 
(red) and black correspond to i = 40° and i = 60° respectively. 



an extended X-ray source is absorbed by a stellar wind, the 
line of sight column density of each X-ray emitting surface 
element is different and the observed flux is given by Eq. [9] 
In general, the observed spectra of a source covered by an 
uniform absorber and of a source covered by stellar wind are 
different (< 5 keV). We assume that the spectral fitting is 
not able to distinguish between different types of absorbers 
if the fluxes at 1 keV are equal i.e. 

1 N 

-^cxp[-a(lkeV)iVH J (j,0)] =expha(lkeV)7Vr x ], (11) 

wher e cr(lkeV) = 2.4 x 1(T 22 cm 2 (|Morrison fc McCammonl 
1983). This is justified, because at energies > 2 keV, the 
photoelectric cross-section decreases quickly and the differ- 
ence between the two absorbed spectra is relatively small. 
Furthermore, the spectrum is most sensitive to absorption at 
energies < 2 keV, but below 1 keV the sensitivity of the in- 
struments starts to drop and the data errors increase which 
may also make it difficult to distinguish between different 
models. 

In order to obtain the lower limit on the size of the 
X-ray emitting region at 1 keV, we explore the parameter 
space (M, 77, Rout, i) to find solutions that satisfy Eg. 1111 
at superior conjunction. In this calculation we reject points 
which are occulted by the star (adjusting N accordingly) to 
test only the effect of absorption. 

The estimated minimum size of the 1 keV X-ray source 
is shown in Fig. [3] The size of the emitting region must be 
larger to compensate for the denser gas when the mass loss 
rate is increased. For 770 — tjmai and at large emitter size, the 
results are not sensitive to inclination changes. The highly 
collimated pulsar wind shock (770 C 1) needs to have a large 
emitting surface to explain the observations. 



4 OCCULTATION IN LS 5039 

The occultation of the X-ray emitting region by the massive 
star does not depend on energy, thus its only signature in the 
spectrum is a periodic reduction in flux. The amplitude and 
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Figure 4. Lightcurves showing the effect of occultation in LS 
5039 (a) i = 60° and 770 = 0.6 and (b) i = 90° and 770 = 0.004 
(measured at periastron). The curves correspond to emitting re- 
gions with sizes R ou t, 15-R* (solid), lOR* (dashed), 5i?* (dotted), 
2R* (dot-dashed). 



duration of occultation depend on three main parameters i, 
rjo and ii ou t. We find that the duration of the occultation 
is longest for binaries with circular orbits and where at all 
phases 770 > 1 (when the CD curves around the massive 
star). The occultation duration can last up to 40% of the 
orbital period. The latter is because curved portions of the 
CD are partly occulted by the star even at phases away from 
the superior conjunction. The amplitude of occultation is 
highest if the inclination is high and the size of the emitting 
region is small compared to the size of the star. However, 
this situation is unlikely to occur in LS 5039 since 770 < 0.6 
(see [Hg. 

When 770 is below unity, occultation influences the shape 
of the X-ray lightcurve only around superior conjunction. 
Lightcurves calculated for the most preferable conditions to 
observe occultation in LS 5039 (high 770, high i and small 
Rout) are shown in Fig.|4^,. The dip in the lightcurve caused 
by occultation is narrow and covers only A(f> ~ 0.2 in phase 
around superior conjunction. The depth of the minimum 
strongly depends on -Rout- Note also that for large R ou t two 
minima appear, separated by a local maximum. This effect 
is related to the size of the shadow cast by the star onto 
the CD surface, which is larger when the line of sight is 
tangential to the surface of CD. These two local minima 
are too narrow (A0 ~ 0.025) to be resolved in the Suzaku 
lightcurve. A similar lightcurve study for 770 <C 1 shows even 
deeper and narrower minimum than in case of 770 ~ 77 max 

(Fig. Hp). 

The X-ray modulation cannot be explained by occulta- 
tion alone but this does not preclude an observable effect in 
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Figure 5. The areas above each curve correspond to permitted 
values of i and i? ut for which effect of occultation at superior con- 
junction does not exceed the 3ad limit. Each curve corresponds 
to different value of 770, solid 0.6, dashed 0.08, dotted 0.02, dot- 
dashed 0.004 measured at periastron. 

the form of a sharp drop in flux around superior conjunction. 
In order to place an upper limit on this effect, we normalize 
the Suzaku lightcurve to 1 and fit it with a sine function in 
order to remove the orbital modulation. The standard de- 
viation of the subtracted data is <7d = 0.1. The upper limit 
on the depth of a dip in the data is assumed to be 3<7d, 
meaning a ^ 30% reduction of the 1-10 keV flux would not 
be statistically significant. The dip in the lightcurve at the 
superior conjunction may be caused by the combined effects 
of absorption and occultation. The effect of absorption on 
1-10 keV flux is low. In a simple estimate, for JVh = ax 
and for the observed power law spectrum of LS 5039 with 
photon index T = 1.51, the integrated absorbed (using uni- 
form absorber, i.e., wabs model in xspec) 1-10 keV flux is at 
most ~ 6% lower than the 1-10 keV unabsorbed flux. There- 
fore we do not take absorption into account in the following 
calculations. 

There is no significant dip in the data of LS 5039 at the 
phases close to superior conjunction. To put constraints on 
the binary parameters based on the lack of occultation only, 
we explore the parameter space to find solution to equation 
in superior conjunction 

1 N 

^X)Ci&0 = l-3o-d. (12) 
j=i 

The final product of the study is Fig[5]where areas above 
each curve correspond to inclination angles and sizes of the 
X-ray emitter for which the effect of occultation at superior 
conjunction does not exceed the 3<7d limit. The range of 
allowed parameters changes with the shape of CD. 



the compact object. Even for the most favorable conditions 
(t}o = Imai ~ 0.6 and large i) the duration of occultation 
never exceeds A<£ = 0.2. The cau se of the X-ray modu- 
lation has to be found elsew here (|Takahashi et al.l 120091 ; 
iDubus, Cerutti fc Henr]|201Ch . 

Constraints on the binary geometry are derived from 
the observed upper limit on the intrinsic column density and 
on the depth of the occultation dip. Limiting the effects of 
absorption and occultation requires the X-ray source to be 
extended. The minimum size of the X-ray emitting region, 
which depends on the stellar mass loss rate, shape of CD and 
inclination angle, varies between 3-15 7?*. An emitter size 
> 4ii* with T]o > 0.02 would be compatible with a 90° binary 
inclination. The limit based on the lack of X-ray eclipses and 
assu ming a point source at the compact object location is 
60° l|Casares et al.ll2005T i. A large X-ray source also loosens 
the constraints on the stellar mass loss rate derived from 
the observed lack of intrinsic absorption. For example, for 
770 ~ 0.004, an emission region ~ 3i?* at inclination ~ 50°, 
allows a mass-loss rate of 1.5 x 10~ 7 Mq yr _1 compared to 
4.7 x 10 -8 Mq yr _1 for a point source located at the pulsar 
position. A more precise determination of the stellar wind 
mass loss rate would greatly help narrow down the possibil- 
ities. Still, we can conclude that an extended X-ray source 
appears necessary to reconcile the pulsar wind scenario with 
the upper limits on X-ray absorption in the wind and the 
absence of occultation features in the lightcurve. 

An extended X-ray source in gamma-ray binary can be 
expected in the pulsar wind scenario. High-energy electrons 
will be accelerated and randomized all along the termina- 
tion shock of the pulsar wind. The shock distance from the 
pulsar is smaller towards the stellar companion than away 
from the O star, the efficiency of particle acceleration may 
also change at different locations. Calculations also show 
that the synchrotron emission from the electrons p eaks in 
the 1-10 keV range only after sig nificant cooling (|Dubusl 
l2006l;lDubus, Cerutti fc Henrill20Qgf ). For a typical magnetic 
field of 0.1-1 G, the injected electrons radiate primarily syn- 
chrotron above 0.1 MeV while they upscatter stellar photons 
to energies above a GeV. The electrons are advected away in 
the shocked flow with an initial speed « c/3. The bulk of the 
1-10 keV radiat ion is emitted at a distance « 3ii* from the 
star in Fig. 4of lDubusl (|2006l l. More detailed modeling is re- 
quired to obtain the exact evolution of the shock conditions 
with distance and quantify precisely the extent of the X-ray 
emission region. Such modeling would also yield the precise 
contribution to the X-ray emission and absorption from the 
shocked stellar wind region, which we have neglected here 



5 CONCLUSIONS 

Using a 3D model of a gamma-ray binary in the framework 
of the pulsar wind scenario, we tested the influence of X-ray 
absorption and occultation on the lightcurve and spectrum 
of LS 5039. We find that occultation cannot be responsi- 
ble for the smooth X-ray orbital modulation in LS 5039, 
since this would require 770 > 1- This appear unlikely be - 
cause there is no sign in the UV lines (|McS warn et al.ll2004T ) 
that the O star wind is quenched on the hemisphere facing 



The application to other known gamma-ray binaries is 
not straightforward since these contain a Be star with a 
dense equatorial outflow in addition to the tenuous stellar 
wind. The geometry of the interaction region with the pulsar 
wind changes between the polar wind and disk wind and can 
have a complicated shape at the transition that has yet to 
be investigated. However, the wider orbits of LS I+61°303 
and PSR B1259-63, the estimated inclination of 36° in PSR 
B1259-63, will limit the impact of absorption and occulta- 
tion on the lightcurve. 
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6 APPENDIX A - RANDOM GENERATION 
OF UNIFORMLY DISTRIBUTED POINTS 
ON THE 3D SURFACE. 

Here we explain how to generate a random sample of points 
uniformly distributed on a 3D surface S, symmetrical with 
respect to x-axis, described by a function /(p) where p — 
^Jy2 _|_ z 2 w e ac j pt an accept-reject algorithm commonly 
used for generating random samples drawn from an arbitrary 
distribution. 

First, we generate N random points distributed uni- 
formly on the yz-plane. The points are denoted in polar coor- 
dinates as (pi,i/>i), (p 2 ,ip2), — , (pn^n) where pi G (0, p max ) 
and tpi G (0, 2ty) . 

Second, we assign to each point a random number 
w G (0, 1), where wi, W2, i»n are drawn from an uniform 
distribution. The Wi are weights, used to accept and reject 
points on the basis of a density function 

P( ) dS = 

[Pi> ~ MpiApAip ~ 



w )' 

where Ap — > 0, -P(pi) G (0, 1), dS is a unit 3D surface ele- 
ment given by Eq. [71 and 



M = ^/(p ma x + Ap)-/( Pmax )^ (M) 

Finally, from our sample of points uniformly distributed 
on the 3D surface, we choose only points for which Wi < 
P(pi). In the paper, / is the solution to Eq. [5] while p ma x 
is determined by Rout and has to be found numerically. An 
example sample of points is plotted in Fig. [2] 



